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■ Non-magnetic impurities affect the paramagnetic response of superconductors via the associated 

' spin-orbit interaction which, when the non-magnetic impurity is ciose to the unitary iimit, must be 

treated beyond the ciassicai Born approximation. Here the Zeeman response of two-dimensional s- 
and d-wave superconductors is calculated within the self-consistent T-matrix formulation for both 
impurity and spin-orbit scatterings. It is shown that at the unitary limit, for which the spin-orbit 
scattering is maximum, the spin-up and spin-down channels becomes decoupled implying full Zeeman 
' splitting of the quasiparticle excitations. These results could be used to test the unitary scattering 

hypothesis in high-T c superconductors. 
PACS number(s): 74.20.Fg, 71.70.Ej, 74.62.Dh 



> 

o 



- I - < 



I. INTRODUCTION 
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The response of a superconductor to defects and/or impurities brings important informations on the nature of the. 
superconducting state and has been the subject of an enormous amount of theoretical and experimental research!!] 
The effect of impurities on the superconducting state depends crucially on the nature of the impurity (magnetic or 

anon-magnetic) and on the symmetry of the order parameter. For isotropic s-wjaizje superconductors, the response to a 
non-magnetic impurity is feeble,cl while the effect of a magnetic one is dramatic.lju Instead when tlie order parameter is 
\ anisotropic, like in d-wave high-T c superconductor, the.response to disorder is always dramatic,ErH leading to Kondo- 
£^ ■ like effects in case of magnetic scattering potentials,!!! or resonant behaviors for non-magnetic impurities close to 
the unitary limit .H For high-T c superconductors therefore it is experimentally more difficult to establish whether the 
impurity acts effectively as a magnetic or a non-magnetic scattering potential. For example, recent scanning tunneling 
microscope images of the local tunneling conductivity around a Zn impurity in Bi2212,B have been fitted by both 
7-H , non-magneticEB and magnetic impurity models J3 

A topic which could be helpful to clarify the effective nature of disorder in high-T c oxides is the analysis of the re- 
sponse to some external applied perturbation. The aim of this paper is to show some important consequences of having 
strong non-magnetic impurities on the Zeeman response of a s- or d-wave superconductor. The way in which non- 
magnetic scattering centers affect th£.scin degrees of freedom is via the associated spin-orbit interaction as described 
. by the so-called Elliott- Yafet theoryOlld Hence, if v is the non-magnetic impurity potentialjlj the corresponding spin- 
y—{ ' orbit scattering is proportional to v so = vSg, where Sg is the shift of the g-factor. The actual value of Sg depends on 
the wave function penetration into the ions and the Fermi surface topology and it is a rather difficult problem.E3 For 
' copper-oxides the main contribution to Sg should come from the d orbital of Cu atoms for which 5g ~ 0.1. Within 
the Born approximation, the spin-orbit scattering rate in the normal state is therefore l/r so ~ (^ff) 2 / T im P <C 1/Ti mp , 
where l/ri mp is the scattering rate due to v alone. This is also known as the Elliott- Yafet relation. However when 
the impurity scattering is close to the unitary limit, as often advocated for impurity doped high-T c superconductors, 
the Elliott- Yafet formula must be generalized in order to include multi-scattering processes. For a two-dimensional 
system with sufficiently jiiluted impurity concentrations rii, the solution of the normal state T-matrix equations for 
(j , both v and v so leads tox3 



JL = 2 1 + C2 — (1) 

l + (2c/^)2 Timp ' { > 

where c = I/ttNqv, l/ri mp = 2r/(l + c 2 ), Y — ih/ttNo and Nq is the density of states per spin direction at the 
Fermi level. In the weak scattering limit c » 1, Eq.(Q) reproduces the result of the Born approximation: l/r so = 
{Sg) 2 /2ri mp <C lri mp . However for c = 0.1, that is the value estimated in Ref. ||, Eq.(|l]) leads to l/r so — 0.4/T imp 
when Sg = 0.1, and in the extreme unitary limit: limc^o 1/tso — 2/T; mp as long as Sg ^ 0. Hence, when the impurity 
potential is strong, or more generally as long as c ^< Sg, inevitably the spin-orbit interaction becomes as important 
as the spin-independent coupling to the impurity. 

The above discussion suggests therefore that if non-magnetic impurities in high-T c superconductors are close to 
the unitary limit, the effect of spin-orbit coupling should be large. In particular, the Zeeman response to an applied 
magnetic field should be deeply altered by the spin-mixing processes associated with v so and eventually, for sufficiently 
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strong spin-orbit scattering, the Zeeman splitting should vanish. Here it is shown that for two-dimensional systems 
this conclusion is actually wrong: the Zeeman splitting resulting from the solution of the T-matrix equation is much 
more robust than that obtained within the Born approximation, and at the unitary limit (c/Sg — 0) both s- and 
d-wave superconductors are fully Zeeman-splitted by an applied in-plane magnetic field. 



II. SPIN-ORBIT T-MATRIX 



For quasi-two dimensional systems the Zeeman response to an external magnetic field H should be best observed 
when H is directed parallel to the conducting plane (for example .the Cu-0 plane in copper-oxides), since in this case the 
coupling of H to the orbital motion of the electrons is minimized.Ej Hence, the total hamiltonian is Ji — Tio+^imp+^so 
where: 

= e ( k ) c L c ka - acl a c ka - Y A ( k )(4T c -ki + c -kjCkT), (2) 

k,a k,a k 

where e(k) is the electron dispersion measured with respect to the chemical potential, a is a spin index, h = p-bH 
and [1b is the Bohr magneton. In the following, it is assumed that the charge carriers are confined to move in the x-y 
plane so that k = (k x ,k y ) and that the spins are directed along and opposite to the direction of the magnetic field, 
fixed to lie along the x-dircction: H = Hx. For s-wave superconductors A(k) = A while for d-wave superconductors 
A(k) = Acos(2(/>) where <f> is the polar angle in the k x — k y plane. Without loss of generality, here A is used as 
an input parameter although it should be calculated self-consistently from a suitable gap equation. Moreover, for 
simplicity, local variations of the order parameter are neglected. The impurity and spin-orbit hamiltonians, 7ii mp and 
H ao , are given by: 

n imp = v £ £ e ~ i(k ~ k,) ' Ric Lc k '«, (3) 

k,k',i a 



= i ^^ e -^- k ')-R i([kxk1 ^ )(c t 

K F .. .., , 



Kso J2 ^ l(k - k,) - R '([k x k'] • z)(c^c kn + 4 ; cv T ), (4) 

k,k' 

where Ri denotes the random positions of the impurities and kp is the Fermi momentum. Note that because of 
two-dimensionality, the spin-orbit matrix element of Eq. (^|) is proportional to a z and, since the spins are quantized 
along the x direction, the spin-orbit scattering is therefore always accompanied by spin-flip processes. The Zeeman 



response for 7Yi mp = and TL so = has already been considered in Ref. |16| for d-wave superconductors and in Ref. 
|l7| for mixed symmetries of the order parameter. The inclusion of 7ii mp , which however does not mix the spin states, 
has been studied in Ref. ^|. The total hamiltonian Tt = H.q + Tti mp + TL so for <i-wave symmetry has been considered 
in Ref. ^ within the Born approximation for both impurity and spin-orbit scatterings. Here, instead, the problem is 
generalized beyond the Born approximation by solving the self-consistent T-matrix equation for both TL lm „ and TL so - 
The generalized Matsubara Green's function G(k, n) in the particle-hole spin space resulting from Eqs.(pH3) satisfies 
the Dyson equation G _1 (k, n) = G^" 1 (k, n) — E(k, n), where G ( ^" 1 (k, n) — iuj n — /Q3e(k) — ,02^2 A(k) — hpsT 3 is the 
propagator resulting from Hq. The Pauli matrices pi and ti (i = 1,2,3) act on the particle-hole and spin subspaces, 
respectively. Whitin the self-consistent T-matrix approach, the self energy is S(k, n) = 7j.;T(k, k, n) where the T- 
matrix is the solution of the following equation: 



T(k, k', n) = «(k, k') + M ( k > k")G(k", n)T(k", k', n), (5) 

k 



and u(k, k') = p^v + i5gv[k x k'] 2 ri. Because of the momentum dependence of the spin-orbit part of u(k, k'), the 
T-matrix can be splitted into the impurity and spin-orbit contributions for both s- and ci-wave symmetries of the 
order parameter. Hence, T(k, k',n) = T imp (n) + T so (k, k',n) where T; mp (n) = p 3 v + p3vJ2]iG(k,n)Ti m p(n) is the 
usual impurity T-matrix and: 

T so (k,k» = i5gv[k x k%n +i6gvJ2$- x k'%r 1 G(k",n)T so (k",k',n) ) (6) 

k" 

is the spin-orbit T-matrix. The solution of Eq.(^|) is of the form0 
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T 80 (k,k',n) = i6gv[kxt(k',n)]r 1 , 



(7) 



where 



t(k,n) = k + iSgvJ2 k ' T iG(k' ,n)[k' x t(k,ra)] 2 



(8) 



The above equation can be easily solved in terms of the components t x (k, n) and t y (k,n) of the vector operator 
t(k, n): 



t x {k,n) =Aj(n) 



t y (k,n) = Aj(n) 



k x + i5gvk y ^2(k x ) 2 TiG{k' , n) 

k' 

k y — idgvk x ^^(fc a ) 2 TiG(k', n) 



(9) 
(10) 



where A x Mn) and A ^(n) are the inverse of the following 4x4 matrices: 



A xy (n) = 1 - {5gvf 



A yx {n) = 1 - (6gv) 2 



^2(k x ) 2 nG(k,n) 

L k 

J2(k) 2 nG(k,n) 



^2(k y ) 2 nG(k,n) 

k 

^(fc x ) 2 n G(k,n) 



(11) 
(12) 



Finally, from Eq.(Q), T so (k, k, n) reduces to: 

T so (k,k,n) = iSgvk x k y [A^(n) - A^(nj\ n + (<5.H 2 £ [^(n)^^) 2 + ^(n)^^) 2 ] r 1 G(k',n)r 1 . (13) 

k' 

Further analysis of the T-matrix problem requires the explicit inclusion of the symmetry of the order parameter. This 
is done in the next subsections where Eq.(n3|) is solved for both s- and d- wave symmetries. 



A. s-wave symmetry 

The usual procedure to evaluate self-consistently the electron propagator is to guess the form of G(k, n) which, after 
being substituted into Ti mp (n) and T so (k, k, n), generates only combinations of pi and Tj matrices already contained 
in G(k, n). The direct substitution of Go(k, n) into T; mp (n) and T so (k, k, n) is a practical way to guess the correct 
form of G(k, n) via the Dyson equation. When this is done, it is easy to realize that when the symmetry of the order 
parameter is s-wave, A(k) = A, the two matrices A xy and A yx defined in Eqs.([ll],[l2|) become equal. Hence, the form 
of the electron propagator for an s-wave symmetry of the order parameter reduces to: 



G 1 (k,n) = i[G) - ihp3T 3 ] - /03[?(k) - iAp 3 r 3 ] - p 2 T 2 [A - if p 3 r 3 ], 



(14) 



where the frequency dependence of the tilded quantities is implicit. The tilded quantities are obtained by substituting 
Eq.(|l4|) into the equations for the impurity and spin-orbit T- matrices and requiring self-consistency via the Dyson 
equation. In general, the solution is very complicated but a considerable simplification arises if infinite electron band- 
width and particle-hole symmetsy of the normal state electron dispersion are assumed. In this case in fact several 
integrals over k average to zeropl leading to the following self-consistent equations: 



iQ± — i(uj n ± ih) + 



1 + c- 



■9± + 2r 



(2c/Sg) 2 g T + g± 



A-i 



1 + c 2 



/± + 2r 



(2c/Sg) 2 f T + f ± 



1 + (2e/6g)* + 2(2c/<S ff ) 2 (/ + /_ - g+gJ) ' 



(15) 
(16) 



where Q± = Cj ± ih, A± = A ± iT and g± = iu)±/ y A± + u) 2 ., f± = A±/y A± + Co 2 .. To display the spin-mixing 
effect of the spin-orbit interaction, equations (|l5 16) are more conveniently rewritten in terms of u± = uj±/A±: 
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OJr, ± ih 



u± 



(1 



-2TTJ2 



(«) +2 (t)7I 



1 + U±Uzfi 



Apart for the trivial limit u± — (uj n ± ih)/A which holds true in the absence of spin-orbit interaction (8g = 0) 
two spin channels u + and w_ are coupled together. Within the Born approximation, c/Sg 3> 1, equation ( J17| ) r 
to the two-dimensional version of the u± formula found in classic literature :EjO 



(17) 



the 
educes 



±ih 



u± 



l T (Sg 
2 \ c 



u± 



(1 



: )l/2 



(18) 



The novel feature displayed by the more general expression ( |l7j ) is that, as the unitary limit c/Sg = is approached, 
u + and U- becomes decoupled and the full Zeeman splitting u + — U- = 2ih/A is recovered. In such a limit therefore, 
the s-wave superconductor is fully Zeeman-splitted as if the spin-orbit scattering would be spin conserving. The same 
conclusion can be obtained by calculating the zero temperature spin susceptibility Xs as inferred by the linear response 
theory. In fact, by including the spin- vertex function consistent with the T-matrix formulation it is possible to show 
that: 



Xs TTT 

Xn A V 1 



1 



1 



K/A)2 [l + K/A) 2 ] 1 /2 + Ps 



(19) 



where % n 

Pso < 1 



Xn = 2/igiVo, T is the temperature and p so — (T/A)(Sg/c) 2 /[l + (Sg/2c) 2 ] 2 . At zero temperature and for 
, Eq.(|l9|) reduces to: 



Xs ^ 

Xn 



7r arccos(p so ) 
2~ v/l-Ps 2 o 



(20) 



For Sg/c <C 1, p S o — (r/A)(<5g/c) 2 and Eq. (|20j) becomes equal to the two-dimensional version of the Abrikosov- 
Gorkov expression based on the Born approximation^ Instead, for c/Sg <C 1, p S o — (r/A)(2c/<5(7) 2 and Xs/Xn — 
2n(T/A)(2c/Sg) 2 < 1 which vanishes when c/Sg = 0. 

The absence of spin-mixing contributions at the unitary limit can be interpreted as a consequence of the fact that 
TCo +7ii mp commutes with S x while TL so commutes with S z . Therefore, for weak spin-orbit scattering {c/Sg 3> 1) S x is 
a rather good quantum number and H. so induces weak spin- flip processes leading to coupled u± equations. The spin- 
decoupling at the unitary limit c/Sg <C 1, could be explained by arguing that for very strong spin-orbit interaction S z 
rather than S x is a good quantum number. The Cooper pairs are then formed by electrons with opposite spins in the 
z direction and the spin rigidity of the superconducting condensate is efficient against spin-flip transitions induced by 
the magnetic field H = Hi. In the limiting case of infinitely strong spin-orbit interaction, therefore H can only induce 
polarization of the quasiparticle excitations. Note that, in case the magnetic field is directed along the z direction, the 
total hamiltonian then commutes with S z and the Zeeman response of a s-wave superconductor becomes independent 
of the spin-orbit interaction for whatever value of c/Sg. Of course, for a three-dimensional system the above reasoning 
does no longer apply because the spin-orbit interaction does not commute with any component of S. 



B. d-wave symmetry 

For the above considerations to be valid it is required only two-dimensionality and a singlet superconducting con- 
densate. Therefore in principle also a two-dimensional <i-wave superconductor should exhibit spin-channels decoupling 
as c/Sg — > 0. This is indeed so even if there are qualitative differences-with respect to s-wave superconductors since 
for d-wave symmetry the spin-orbit scattering becomes pair breaking.E3 Again assuming particle-hole symmetry and 
an infinite electron band-width, for A(k) = A(<fi) = Acos(2</>) the electron propagator is of the form: 

G^iKn) = i(G> - ihp 3 T 3 ) - p 3 [e(k) - iAp 3 r 3 ] - p 2 r 2 [A(0) - if(0)p 3 r 3 + inQ^)], (21) 

where A(<fi) = Acos(2</>), T((f>) = Tcos(2(p), and Cl(<p) — f2sin(20). The origin of (l(4>) (absent in the s-wave case) 
stems from the fact that, for e?-wave symmetry, the two matrices A xy and A yx in Eqs . (JT^JT^) are no longer equal, 
so that the term proportional to k x k y = sin(2</>) in T so (k, k, n), Eq. (|l3|) , is nonzero. As for the s-wave case the 
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self-consistent Dyson equation can be expressed in terms of Q± and A±, but now there is an additional equation for 
( l.'U j r o , or (2c/Sg) 2 g T + (f^-g 2 T )g ± 

luJ± = zK ± lh) + c^(gW 9± + ™[ww- 9+9 --f + f-Y-i 9+ f-+9-ur (22) 

A a , or (2c/,5 g ) 2 / T -(/|~g 2 )/ ± 

_ [{2c/5g) 2 _ _ /+/ _]2 _ Cg+/ _ + g_ f+) 2 > 

q op (2c/<5ff)(ff+/- +g-f+) . . 

[(2c/^) 2 -9 + ff--/ + /-] 2 -(.9 + /-+.g-/ + ) 2 ' 1 j 



.9± 



2 y^ sin(0) 2 ig±jft| + -Eip(k) 2 ] T 2j(g+ - a_)f2(^) 2 

7TiV ^ 2^ [itf +£+(k) 2 ][c£ +^_(k) 2 ]-^(,/)) 2 {(^ + -^) 2 + [A + (0)-A_(^] 2 }' 



, 2 v sin(0) 2 A ± (0)p 2 + ff T (k) 2 ] T [A+W - A_(0)W) 2 

J± TT^o^ 2^ +jB+(k )2]p2 + ^_(k)2]-^(0)2{(^ + _^_)2 + [A + (0)_A_(0)] 2 }' 



where £±(k) 2 = e(k) 2 + A±(</>) 2 + fi(0) 2 and g° ± can be obtained from Eq.Q by setting sin(0) 2 -» 1/2. The off- 
diagonal contribution f2 defined in Eq.(p4|) is responsible for spin-mixing terms appearing in g±, f±, and . However, 
at the unitary limit c/Sg = 0, vanishes and the above self-consistent equations reduce to: 

iui± = i(u n ± ih) - -Q- + 2r 2 g± 2 , (27) 



where 



2 i# fq) ± BinW 2 

' 2 7 r[A ± (0) 2 +w 2 : ] 1 /2' v ; 



. # A ± (0)sin(0) 2 
/± = 2/ ^[A ± (0) 2 + ^ 2 ]V2- ( 3 °) 



The two spin channels + and — in Eqs.( |27| , |28| ) are now completely decoupled in analogy therefore with the s-wave 
case treated before. However now even at T = the spin susceptibility is expected to remain non-zero (as long as 
r / 0). This is due to the fact that, although there are not spin- mixing processes at c/Sg = 0, the pair- breaking effect 
of both impurity and spin-orbit scatterings leads to a finite density of states at the Fermi level. Ej In this situation 
therefore the Zeeman splitted density of states is a more direct evidence for the spin decoupling effect at c/Sg — * 0. 
This is shown in Fig. 1 where the two spin channels density of states, N±(u>), are plotted for T = 0.1, Sg = 0.1, 
h = 0.2 and for different values of c. N±(u>) is calculated numerically from: 

^ = -sgnMlm[ ff °M], (31) 

where g±(to) is the analytic continuation on the real axis (ioj n — > oj + id) of Eq.(p5|) (with sm(c/)) 2 — > 1/2). @ In Fig. 
la, N±(uj) is calculated from the solution of the general equations (|22|)-(|2~i|), while, fop-comparison, the result for the 
Born approximation applied to the spin-orbit part of Eqs.( p2[ - |24| ) is shown in Fig. lbE3 Up to c = 0.1 (5g/c = 1) the 
general T- matrix solution and the Born approximation agree quite well, while already for c = 0.05 {Sg/c — 2) the 
splitted coherence peaks at cu/A ~ ±(1 ± h) are still quite visible in Fig. la and completely suppressed in Fig. lb. 
Since at the unitary limit c/Sg = the spins are completely decoupled, Eqs.( |27| , p8| ), the two spin density of states 
are identical and shifted by ±h one respect to the other. This is in contrast to the Born solution, Fig. lb, for which 
the strong spin-mixing terms lead to a flat density of states. It should be stressed that before reaching the Sg/c^> 1 
limit, the Born approximation predicts that d-wave superconductivity is already completely destroyed.! 2 !! Therefore 
the results for c = in Fig. lb should be considered just as a mathematical limit to be compared with the solutions 
of the T-matrix approach of Fig. la. 



5 



III. CONCLUSIONS 



In summary, it has been shown that when non-magnetic impurity scattering is close to the unitary limit, the associ- 
ated spin-orbit interaction is not small provided c ^< Sg and must be treated beyond the simple Born approximation. 
Within the self-consistent T-matrix approach, it has been demonstrated that the Zeeman splitting of both s- and 
ci-wave two-dimensional superconductors is much more robust than that obtained by the Born approximation. At the 
unitary limit c/5g = 0, for which the spin-orbit coupling is maximum, a two-dimensional s-wave superconductor does 
not show any spin-mixing processes and the Zeeman response coincides with that of a pure superconductor. Also for 
d-wave superconductors the spin-orbit coupling becomes effectively spin-conserving at c/<5e-j= 0, but in addition it 
induces pair breaking effects which must be added to those caused by the scalar impurities JIj 

Let us comment now on the possible limitations of the present theory. The calculation method here used is a standard 
one based on a T-matrix approximation for diluted impurities JjBES However when applied to two-dimensional systems 
like the-eopper-oxides, the standard procedure is complicated by the appearance of singularities in the electron self- 
energy.cSI Contraxy to the results based on the T-matrix solution and on self-consistent approaches to deal with 
the singularities]^ non-perturbative methods suggest that the density of states of a d-wave superconductor actually 
vanishes non analytically at the Fermi level.E3 The low energy behavior appears to be heavily modified by the 
level spacing, of a localization volume which leads to the opening of a pseudogap in the low- lying single electron 
excitations.Eil It should be however noted that discepancies between different approaches affect only the very low- 
energy excitations, while for energies not much smaller than A the T-matrix approach is quite realiable (for finite 
but small impurity concentrations). In this respect, the main result of Fig. la (i. e. the persistence of the Zeeman 
splitting of the coherence peaks in the density of states even when c/Sg is zero) should not be an artificial feature 
of the T-matrix approximation. This conclusion is also sustaine d by the quite general physical explanation of the 
Zeeman-splitting persistence at the unitary limit proposed in Sec. II A and from an analysis of the singee spin-orbit 
impurity problem not reported here. Note that, for these same reasons, some standard simplifications employed in the 
present calculations (infinite band- width, particle-hole symmetric electron dispersion and absence of local suppressions 
of the order parameter) should not affect too seriously the main result. 
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FIG. 1. The Zeeman-split quasiparticle density of states N+(u))/No (dashed lines) and N-(lo)/No (solid lines) for a d-wave 
superconductor with h = H/A = 0.2, F — 0.1, <5g = 0.1 and different values of the scattering parameter c. (a): solution of the 
complete T-matrix equations, (b): solution for the Born approximation to the spin-orbit coupling. 



8 



